Surface modes of ultra-cold atomic clouds with very large number of vortices 
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We study the surface modes of some of the vortex hquids recently found by means of exact 
diagonahzations in systems of rapidly rotating bosons. In contrast to the surface modes of Bose 
condensates, we find that the surface waves have a frequency linear in the excitation angular mo- 
mentum, M > 0. Furthermore, in analogy with the edge waves of electronic quantum Hall states, 
these excitations are chiral, that is, they can be excited only for values of / that increase the total 
^ angular momentum of the vortex liquid. However, differently from the quantum Hall phenomena 

■^1^ , for electrons, we also find other excitations that are approximately degenerate in the laboratory 

frame with the surface modes, and which decrease the total angular momentum by I quanta. The 
surface modes of the Laughlin, as well as other scalar and vector boson states are analyzed, and their 
observable properties characterized. We argue that measurement of the response of a vortex liquid 
to a weak time-dependent potential that imparts angular momentum to the system should provide 
valuable information to characterize the vortex liquid. In particular, the intensity of the signal of 
the surface waves in the dynamic structure factor has been studied and found to depend on the type 
, of vortex liquid. We point out that the existence of surface modes has observable consequences on 

ly-j ' the density profile of the Laughlin state. These features are due to the strongly correlated behavior 

' of atoms in the vortex liquids. We point out that these correlations should be responsible for a 

remarkable stability of some vortex liquids with respect to three-body losses. 



PACS numbers: 3.75.Fi, 05.30.Jp, 73.43.-f 



I. INTRODUCTION 



O ' Everybody is familiar with the phenomenon of vortices. We observe it every time we watch water down a drain. 
In classical fluids, it is a consequence of broken Galilean symmetry 1']: A fluid that fills a region of space provides 
us with a privileged reference frame, namely the one where the fluid is at rest. Nature tries to restore the full 
. Galilean invariance at the vortex core, which explains why the fluid density drops to zero (or almost zero) there. 
^ ' In superfluids, however, vortices exhibit certain peculiarities not seen in classical fluids. This is because quantum 
coherence is maintained throughout the superfluid volume, and since the wave function must be singly valued the 
circulation around a vortex is quantized. Therefore, the stability of vortices is ensured by topological rather than 
dynamical reasons Nevertheless, a vortex carrying n > 1 circulation quanta is unstable with respect to decay into 
' n singly quantized vortices. When many of these have appeared, they form a triangular lattice, the Abrikosov lattice. 
fSI In ultra-cold atomic gases the formation of this lattice has been observed in a recent experiment Q . If the rotation 
. frequency, fi, is further increased to approach the trap frequency Q, several authors [ESHl have pointed to the 
interesting possibility that the Abrikosov lattice melts due to quantum fluctuations. This regime {fl ^ is known 
as the critical rotation limit, and some experiments have already begun to be explore it 0|. 
y ' What kind of phenomena can be expected to emerge in the critical rotation limit? Under certain conditions, 
' , I which will be discussed below, it has been predicted that bosons organize themselves in highly correlated, two- 
^ • dimensional, "vortex liquids". These are states that cannot be described by the standard (i.e. Gross-Pitaevskii) 
Q mean field theory 5, 10] . Instead, they seem to be quite accurately described by microscopic wave functions la, llOl , 
Q closely related to those used for electrons in the context of the fractional quantum Hall effect (FQHE) [Til Il5j |. 
^ ' Recent numerical studies using periodic boundary conditions |^ have shown that, just like their FQHE counterparts, 
• i-H , homogeneous vortex liquids are incompressible. This means that changing their density requires a finite amount 
of energy, signaling the existence of a spectral gap. When carried out for small systems in a harmonic trap, exact 
H , diagonalization studies |0, 1^ found a series particularly stable states at "magic values" of the angular momentum, 
- - I some of them related to the homogeneous vortex liquids Q. 

In this work we address the following question: Which are the experimental signatures of these vortex liquids? 
The question is relevant since experiments driving ultra-cold atom clouds to the critical rotation limit are likely to 
proliferate in the near future. It is also important to find ways to characterize the vortex liquids by non- destructive 
means, especially if one intends to use their entanglement properties for quantum computing purposes T^. Below 
we shall try to answer the previous question by studying the properties of surface excitations of some of these states. 
Since the vortex liquids are effectively two dimensional, by "surface" we mean the one- dimensional boundary in the 
plane perpendicular to the rotation axis. Excitations along the rotation axis will occur at much higher energies and 
are therefore frozen at the low temperatures of interest here. As we shall argue in section IhTI deformations of the 
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FIG. 1: Rapidly rotating cloud of ultra-cold atoms in the critical rotation regime Q ^ lu±. 



boundary correspond to the lowest energy excitations of the system. Their spectroscopic analysis should provide 
valuable information on what kind of vortex liquid one is dealing with. In the case of the FQHE exhibited by two- 
dimensional electron gases under strong magnetic fields, characterization of the state of the sample is usually achieved 
by means of transport measurements. In particular, a measurement of the Hall conductance automatically yields the 
filling fraction v, which in the present context and in the limit of large particle and vortex numbers can be defined 
as the ratio of the number of particles to the number of vortices, i.e. v = N/N^. However, nothing like transport 
spectroscopy is yet available for ultra-cold atomic clouds. Instead, experiments where the trapping potential is weakly 
deformed in a time-dependent fashion, therefore exciting surface modes, seem more feasible. That this is feasible with 
a rotating cloud containing a large number of vortices has been already demonstrated in a recent experiment 
(however, not yet in the regime in which we are interested here). 

The study of the surface modes is also interesting because, at low temperatures, these excitations should dominate 
the properties of the system that can be probed by weakly coupling to it. Furthermore, it has been shown in the 
context of the FQHE 19, that analyzing the boundary excitations is the way to characterize the bulk states (i.e. 
their quantum orders 20]). What is more, the boundary of a vortex liquid is also a clean example of an exotic type 
of quantum liquid: The chiral Luttinger liquid. The main difference with the electronic FQHE is that the constituent 
particles of the system are neutral bosonic atoms, which interact via short-range interactions. This setup offers some 
advantages over the two-dimensional electron gas at high magnetic field. The systems are intrinsically clean, and the 
absence of long-range interactions does not lead to some complications introduced by the Coulomb potential (e.g. 
edge reconstructions 'l5] ) . One can also consider atoms with internal degrees of freedom, which opens the possibility 
to study novel quantum Hall states 0,|23. The disadvantages are the fragility of the vortex liquids, which will 
require careful experiments with small numbers of particles, iV < 100 (although many replicas of a rotating cloud can 
be created in an optical lattic e 1211 ). In the future, other methods to "simulate" high magnetic fields using optical 
lattices may become available |22j . and this could lift the constraint on the particle number. However, large systems 
confined by smooth potentials, like the ones used to trap ultra-cold atoms, should exhibit more complicated boundary 
structures like "composite edges" with many branches of surface modes. Nevertheless, the theory developed here 
should provide the basis for the understanding of such systems. 

The experimental achievement of the vortex liquid states to be discussed below is challenging, but on the theoretical 
side a complete understandin g of their properties also poses many challenges. The existing work has mainly focused 
on ground states of scalar |lOll and higher-spin bosons 0, |^ , but the situation is far from being as clear as 
in the FQHE for electrons T5]. In this paper, we shall consider the surface excitations of several ground states of 
scalar and vector bosons. One of them, the Laughlin state, is an exact ground state for rotating scalar bosons. The 
other ground states for scalar bosons that we will consider are approximate, and have been found to exhibit good 
overlap with the exact ground states obtained by exact diagonalization methods 0, 0, ^| . We shall consider the 
Laughlin state first, for which the microscopic construction of the surface waves will be presented. An effective field 
theory will be subsequently developed and shown to agree with the microscopic theory. Althou gh t he assumptions 
for the effective theory strictly hold in the large N limit, it is well established numerically (e.g. [l9l 133. Is^l that it 
also applies to small systems with N ~ 10, within the experimental reach in the near future. 
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As many of the concepts and methods introduced and used in this paper should be new to the community working 
in ultra-cold atom systems, we have adopted a pedagogical approach at the cost of producing a longer paper. We have 
also tried to make the article as self-contained as possible, but without omitting references to the original literature on 
the quantum Hall effect, which should be consulted for more extended explanations. The paper is organized as follows: 
In the next section, the experimental conditions for the existence of vortex liquids are discussed. Next, the simplest 
(from the theoretical point of view) type of vortex liquid, namely the Laughlin liquid, is considered. In section IIVI 
surface excitations of more complicated (but approximate) states of scalar bosons are studied. In the following section, 
we take up vector bosons and analyze the surface excitations of their singlet states. Finally, the conclusions of this 
work can be found in Sect. I VII In the appendices we include some supplementary material. Thus appendixIXldiscusses 
the relation between the laboratory and rotating frames, and how measurements in the laboratory frame are related 
to calculations in the rotating frame. We also discuss the generalized Kohn's theorem in this appendix, and how it 
constrains the energy and peak intensity of the dipole mode. Appendix^ however, gives a brief introduction to the 
plasma formalism for quantum Hall wave functions. More importantly, we also show there how to obtain the wave 
functions from correlation functions of the field theory that describes the boundary excitations. 



II. EXPERIMENTAL CONDITIONS FOR THE OBSERVATION OF VORTEX LIQUIDS 

Before proceeding any further, it will be useful to recall the conditions under which it is expected that a rotating 
cloud of ultra-cold atoms will become a vortex liquid. The strong analogy between the quantum mechanics of rotating 
particles in a harmonic well and charged particles moving in two-dimensions under a strong magnetic field will be 
useful here; therefore, the appearance of Landau levels is expected. To see this, notice that the one-body part of the 
Hamiltonian in the rotating reference frame (see Appendix ^ , 
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can be written in two equivalent ways: 



,1) _ (p,-Afflzxr,)^ ^ M. , , ^ pI, ^ M^ii , 



and 



i^f ) = + (., - f^) z . L, + + ^z^. (3) 



We have assumed an axially symmetric trap with the rotation axis coinciding with the symmetry axis of the trap; the 
following notation has been introduced: Vj = {xj , yj ) and pj — (p^j , Pyj ) are the j-th particle position and momentum, 
respectively, while Lj — Vj x pj is the angular momentum along the rotation axis. In both forms of Hj, it is manifest 
that particles, in the rotating reference frame, feel an effective magnetic field directed along z, whose magnitude is 
proportional to either MQ or Mu;±, depending on the form, M being the atom mass. Perhaps the first form, H^^\ 
shows more explicitly what happens as fl approaches uj±, and the system enters the critical rotation regime: First, 
notice that particles under rotation experience a reduced confinement in the plane perpendicular to the rotation axis. 
This means that, as fl is increased towards the trap frequency more and more particles will accommodate into 
the lowest axial level (this can be helped by making W|| > ujj_ [J). Eventually, when all the particles are in this level, 
the system effectively becomes two-dimensional. At the same time, increasing increases the strength of the effective 
magnetic field and therefore, for sufficiently weak interactions, all particles will make its way into the lowest Landau 
level (LLL). 

To study the quasi two-dimensional system in the critical rotation regime with all the particles lying in the lowest 

(2) 

Landau and axial levels, it is more convenient to use the second form of the Hamiltonian, iJj . Dropping the axial 
part and introducing the following operators: 



(4) 

TTj = '^jx ^'^jyt (^) 

where iTjx = Pjx + Mcj^yj and Wjy — pjy — Muj_\_Xj^ and [tt^, tt^] = AMTiuj_\^^ allows us to write: 

Hf^ = ^ + i^^-^)L3+^^- (6) 
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FIG. 2: For non-interacting particles, the large degeneracy of the lowest Landau level is lifted only by the confining potential. 
In the reference frame that rotates with the system, a "band" of width « TivN/v is formed. 



One can easily convince oneself that this form is diagonal after noticing that ttj and 7fj are nothing but ladder 
operators, which move particles between consecutive Landau levels: tTj promotes the j-th particle to a higher Landau 
level whereas tt^ undoes this operation. Thus, if a particle is in the LLL, the corresponding -Kj must annihilate the 
state. Once diagonalized, one finds that the single-particle orbitals in the LLL have energies (in the rotating frame) 
equal to e(l) = huj± + hvl, where v = — f2) <C uj± and M is the angular momentum. They form a "band" (see 
Fig. I2J) within the LLL. The orbitals themselves take the form: 

^,(.,y) = -l=f^Ve-WV2^^ (7) 



where i = y^h/M u}± is the oscillator length of the trap, and z = x + iy denotes the position of the particle on the 
plane perpendicular to the rotation axis. As the system is effectively two dimensional, we have omitted the orbital 
describing the motion perpendicular to the plane. 

We are now ready to take up the discussion of the conditions under which the scenario described above can 
be experimentally realized. Interactions between the particles are needed for the stability of the vortex liquids to 
be discussed below. Taking them into account, the total Hamiltonian reads H = X^^Li + Si<j=i ^iii ^ij — 
g S'^^'^ {vi — rj)S{zi — Zj), where g is related to the s-wave scattering length, a, as usual: g — Anh^a/M. The interactions 
should be neither too weak nor too strong. Since hv is the level spacing of the single-particle band in the LLL, we 
must have that gn ^ hv but gn < Nhv/v, where n is the mean density of the system. The last condition is the 
requirement that the width of single-particle band must be larger than typical interaction energy, gn. The filling of 
the band is determined by the total angular momentum through the filling fraction v « N"^ /2L. Hence N/v « 2L/N 
measures the mean angular momentum per particle. Finally, since all the particles must be in the LLL, mixing with 
higher Landau levels is avoided provided that Nhv/v <C UQ, and gn <C Kl and at temperatures T ^ h^. Although 
these conditions look very demanding, we expect the rapid progress that has characterized the field over the last years 
will make vortex liquids experimentally available in a near future. 



III. SURFACE MODES OF THE BOSONIC LAUGHLIN STATE 



We begin by considering the simplest case, from the theoretical point of view. A droplet of N bosons in the Laughlin 
state is described by the following wave function 

N 

$^(zi,...,ZAr)= n {Z,-Z,r, (8) 
i<3 = l 

where m = 2 (m must be even for bosons and odd for fermions). In this expression we have only written the polynomial 
part of the wave function. We have omitted (and will omit henceforth) a factor that includes the normalization 
constant and a product of single-particle orbitals involving the gaussians that depend on \zi\^ and the orbital that 
describes the state along the rotation axis. By demanding all particles to be in the LLL, the above function can only 
depend on Zj, and not on its complex conjugate, Zj = Xj — ii/j. To understand this physically, notice that in the 
LLL the kinetic energy is minimized while the angular momentum is maximized. Therefore, the angular momentum 
of every particle along the rotation axis must be positive or zero. This is true as long as the polynomial part of 
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the wave function depends only on cx e^'^J , with / > 0. In other words, the polynomial part of the wave function 

must be an analytic function of Zj, for j = 1, . . . , N. Therefore, allowing to act on <i>„i yields N times the 

energy of a particle in the lowest Landau and axial levels, h(uj± + a;||/2), plus the confinement energy, vL^, where 
Lo = fimN{N — l)/2 = hN{N — 1) is the total angular momentum of the Laughlin state (to see this, consider a 
rigid rotation of the system where Zi Zic"^^ , the angular momentum can be read off the phase factor of the many- 
body wave function). Furthermore, one can easily see that the Jastrow structure of ^ produces a state with zero 
interaction energy for Vij = om S'--^"^ (r^ — fj), which in this case represents the relevant interaction for ultra-cold atoms 
confined to two dimensions j4Q| . 

Let us now consider the excitations about the Laughlin state. Bulk excitations 0| come in two types: quasi-holes 
and quasi-particles. They correspond to a deficit (quasi-holes) or an excess (quasi-particles) of ^ = ^ boson in the 
bulk of the droplet. Another interesting property, which will be useful in the following section, is that creating a 
quasi-hole increases the average angular momentum of the system, whereas creating a quasi-particle decreases it. For 
example, by creating a quasi-hole at the center of the a droplet: YiiLi ^i^m, so that the resulting state is an eigenstate 
of L, the angular momentum increases by hN , and the energy (in the rotating frame) by fivN . However, if we consider 
wave functions of the form [23.12^ s/$,„, where the factor is a symmetric polynomial, s/ = X^iLi with < Z <C iV, 
the angular momentum increases only by hi and the energy by Tivl, much less than the quasi-hole energy. Thus we 
see that the states created by multiplying the Laughlin state by (products of) s; are indeed the low-lying excitations 
of the system. These states describe deformations of the droplet boundary, which in the ground state has a circular 
shape. This is true for all / > 1 except for Z = 1, as si = X^i^i represents a small translation of the droplet center 
of mass, which leaves the shape of the boundary unchanged. 

On the other hand, quasi-particles are not the only kind of excitations that decrease the angular momentum of 
the system. Consider the states created by the action of the operators di = Tl!j=i (^ = 1j 2, . . .) on the Laughlin 
state, where Ttj has been defined in Eq. jSJ. According to the discussion that we made in previous section, it is easy 
to understand what these operators do: They move up one particle by I Landau levels. At the same time, one can 
show that they decrease the total angular momenta by I quanta. The one-body part of the Hamiltonian thus yields an 
excitation energy (in the rotating frame) equal to h(LU± + Q)l, but strictly speaking the states are not exact eigenstates 
of the full many-body Hamiltonian since they do not diagonalize the interaction potential. However, their interaction 
energy for I ^ N is much smaller than their one-body energy since gn <^ hfl, and therefore they can be considered 
as good approximations to the exact eigenstates. Furthermore, one can show that for I = 1, di ^ 'Ylfj=i^j ~ ^i' 
which has zero interaction energy. Thus di and si are the two independent modes that describe the center-of-mass 
motion in the plane perpendicular to the rotation axis. They are sometimes known as "Kohn modes", and their 
properties are discussed in detail in appendix fXl Returning to the general properties of the states created by c?;, and 
applying the relation Ehot = E^ab — (see appendix we find that the energy of these states in the laboratory 
frame is « huj±l, which means that they are practically degenerate (in the laboratory frame) with the surface modes 
described by s;. This forces us to clarify what we meant above by "low- lying" excitations of the system. First of 
all, it is necessary to recall that the system is at thermal equilibrium in the rotating frame and not in the laboratory 
frame. The characterization of excited states as low-lying or otherwise high-energy states must be done in this frame, 
where the states di^m have much higher energy than the states si^m- Therefore, at the temperatures T <C hfl of 
interest here, the statistical probability of finding the system in the latter states is much higher than in the former, 
and therefore the low-temperature properties are dominated by the states generated by s/. Nevertheless, this does 
not mean that if the system is probed by an external field of frequency uj « luj_1 [l — 1,2, . . .) which excites modes of 
either chirality with equal probability, these excitations are not important to understand the response of the vortex 
liquid. As the frequency of the external field is determined by an apparatus located in the laboratory frame, modes 
of both chiralities will be excited. This is a distinct feature of vortex liquids compared to the physics exhibited by 
electronic quantum Hall systems. However, as it will be shown below, the response to an external probe is weighted 
by a thermal factor [1 — exp (— ?ia;flOT(0/^)]~^j ^iid at temperatures such that hfl T fiv, which should be 
more easily accessible, it yields a much higher intensity for the surface modes s/ than for the di modes. Actually, the 
considerations just made regarding the di modes also apply to more general (vortex-liquid) states than the Laughlin 
state. With these remarks, we close the discussion of the excitations generated by the operators di, and we shall not 
discuss them any further in this work. Instead, we shall focus on the surface modes generated by s/ , which exist only 
when the total angular momentum is increased by a small number of quanta I >Q, and for which all particles remain 
in the LLL. 

In order to compute the low-temperature properties of a droplet of vortex liquid, it is convenient to develop an 
effective field theory that describes the low-energy part of the spectrum. The states generated by the action of s/ on 
$m have a number of important properties which must be captured by the effective theory. First of all, they are all 
chiral since I > and cannot be negative for otherwise the resulting state would not be in the LLL. Second, the linear 
dependence on I of the energy implies that the spectral degeneracy is given by p{l), the number of partitions of I. This 
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is the number of distinct ways I can be written as a sum of non-negative integers. To understand this, recall that the 
angular momentum, and hence the energy, are proportional to the degree of the homogeneity of the polynomial part 
of the wave function. Therefore, given that si^m has an excitation energy equal to hvl, so do the states s;_isi$m 
{I = {I — 1) + 1), or s;_2(si)^*&m {I — — 2) + 1 + 1), etc. In the following, we will show that these features are in 
fact captured by quantizing a two-dimensional liquid drop model [l9.] . 




R r 



FIG. 3: (a) Deformation of the droplet boundary and (b) semiclassical approximation to the ground state density profile. 

Following Wen [19| . we develop an effective low-energy description for a finite droplet of = ^-Laughlin liquid. 
Henceforth, we shall consider only the Hamiltonian in the rotating frame. However, when computing physical prop- 
erties we shall transform to the laboratory frame (see Appendix for details of how this transformation must be 
performed). In the ground state, the shape of the droplet is circular, and the density will be taken to be uniform 
and equal to po = N/A, A — ttE? being the area of the droplet (see Fig. ISJ. We find A by noticing that if $,„ is 
expanded in powers of, say, zi, the highest power is /max — fn(N — 1), which corresponds to the angular momentum 
of the highest occupied single-particle orbital in the LLL. The extent of this orbital gives the radius of the droplet, 
R = 'max i ~ (mN)^^^ £; hence po = l/{'rm:£'^) for iV 3> 1. Next, imagine a deformation of the droplet such that the 
boundary is shifted to i? -I- h{6) for < 6 < 2tt. Clearly, since this is a low-energy description the deformation must 
be small: h{9) <C R- In addition, our hydrodynamic description necessarily breaks down at the scale of the oscillator 
length i. This implies that h{9) ^ i or, in terms of angular momentum of the excitations, the effective theory will be 
valid for modes with angular momentum ^ hyrnN . 

In what follows, we work with the particle density per unit angle (sometimes we shall call it current because it is 
proportional to it), 

j{e) - poRh[e). (9) 

Let us write down the energy cost (i.e. the Hamiltonian) for a deformation of the droplet, being v{y) = M{lo\ — 
VP') r^/2, the confining potential in the rotating reference frame (cf. Eq. ||2Jl), and p, = u(|r| = R), the chemical 
potential. The appropriate expressions read: 

H = jd\ {v{v) - p) [p{v) po(r)] - ^ fj de h^io) = ^ £^ defie), (lo) 

where Po{r) = PoO{R — r) and F = v'{R). We will allow for the possibility that the deformation involves the addition 
or removal of a small number Q = f^^ d9 j{9) <^ N oi particles at the boundary. Canonical quantization of can 
be achieved with the help of the continuity equation, 

dtj{l,t)^zvlj{l,t), (11) 

which has been expressed in terms of the Fourier modes of j{9) = X); j(0 "^^P {—US) /2tt. Notice that, contrary to the 
case of single-particle orbitals, in this case / measures, not the absolute angular momentum of the excitation, but the 
increase in the angular momentum over the ground state. The quantized effective Hamiltonian takes the form 

The canonical commutation relations lead to the following current algebra: = -^Si+i'fi, which in the 

mathematical Physics literature is known as U(l) Kac-Moody (KM) algebra. Nevertheless, both the Hamiltonian 
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and the current algebra can be written in a more familiar way if, for I > 0, we introduce b{l) = —iy^Ji—l) and 
b^{l) = i-\/^j(0; hence [b{l),b^ {I')] = 6i^i', i.e. the familiar "phonon" algebra. In terms of these operators 

H = Y,^vl bHl)b{l) + Y^Q^^ (13) 

i.e. the Hamiltonian of a system of chiral phonons with energy dispersion hvl {huj±l in the laboratory frame). 
Therefore, the spectral degeneracies are given by p{l). The chirality of surface modes stems from the broken time- 
inversion symmetry existing in this situation: The system rotates in one direction, which is reversed under time 
inversion. Furthermore, it is also important to point out that the state annihilated by the b{l) operators must be the 
Laughlin state. This is related to the existence, in the rotating frame, of a finite energy gap which must be overcome 
in order to decrease the angular momentum below that of the Laughlin state, and should be interpreted as due to 
the incompressibility of the Laughlin state (i.e. the droplet cannot further shrink) '33]. Thus we conclude that the 
effective theory correctly describes all the features of the low-energy spectrum obtained using wave functions. 

The first property that we want to compute using the effective theory is the response to a weak time-dependent 
deformation of the the trapping potential. This couples to the particle density at the boundary, and therefore to j{0) 
(see appendix 1X1 for more details). The measurable response will be given in terms of the dynamic structure factor. 



+C30 /'27r 

iW 



{l,uj)= dte'^' d0e-*"'(j(e,i)j-(O,O)), (14) 



where the brackets stand for thermal average over the canonical ensemble in the rotating frame, which is where the 
system is in thermal equilibrium. The above function can be written in terms of the response to a external potential, 
by virtue of the fluctuation-dissipation theorem |39| |. At a temperature T, the relationship is given by 



2nimx(^t^) 



where xi^^) is the Fourier transform of 



x(0,i) = -J^W([j(^,i),j(O,O)]), (16) 

where "dit) is the step function. The previous correlation function can be readily computed using current algebra 
described above and the continuity equation, which implies that j{0,t) ~ j{9 — vt). Thus, applying the results of 
appendix we arrive at the following result for the dynamic structure factor in laboratory frame: 

Slab{1.^) ^ Sb.ot{1.^- l^)^^ ^ _^l-n(Z-n)i/T - (1^) 

Since there is a single branch of phonons, the fact that S'lab(^,'-^) is peaked at a; = ujj_l should not be surprising (in 
real systems, a finite broadening of the peak is expected). However, it is also important to notice that the spectral 
weight of the peak associated with the surface mode (cf. Eq. H17|) ') is proportional to the filling fraction i/ = 1/m oi 
the Laughhn state. 

A way to understand the proportionality to the filling fraction of the spectral weight is to show that, as far as the 
response to an external potential is concerned, a droplet of Laughlin liquid behaves as a rapidly rotating cloud of 
non-interacting particles (i.e. an ideal gas) with peculiar exclusion statistics p^ . To understand what this means, 
let us assume that when the ideal gas is in thermal equilibrium it is described by a density matrix p° such that 
p^\k) = N{k)\k). The function N{k) is the occupancy of a single-particle orbital |fc) with angular momentum equal 
to hk. To compute the linear response of this ideal gas to the external perturbation Svext{0,t) localized near the 
boundary, let us consider the linearized equation of motion for the perturbed density matrix -|- 5p{t), 

th{k + l\dtSpit)\k) = {k + l\ [Ho, 6p{t)] + [Sv^xtiO, t), p°] |fc). (18) 

where Ho\k) — e{k)\k), e{k) = hvk being the single-particle dispersion. We thus find (r/ 0+): 

n ^= V ^WM) _ 1 V iV(fc)-iV(fc + 

Sv'I^Lo 5vea^t{l,uj) 2nh ^uj-{e{k + l)-e{k))/h + i7]' ^ ^ 
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Since the particles have been assumed to be non-interacting we can use the results of Sect. which imply that the 
single-particle dispersion is linear with the angular momentum. Thus the previous expression reduces to 

= — Y[N{k)- N{k + l)]^ ^ . (20) 

' ^ 2TThuj-vl + iri f^^^ ^ ' ^ 2TTh uj - vl + iri ^ ' 



It is important to stress that the above expressions make sense as long as I ^ \/rriN, otherwise the effect of the 
external potential is not reduced to the neighborhood of the boundary and radial part of the single-particle orbitals 
must be taken into account. In the previous expression we have denoted as g"^ the mean occupancy at the bottom 
of the single-particle band (cf. Fig. |2Jl, which at the low temperatures where the gas is quantum degenerate must be 
a constant. By comparing this expression for xih^) with the one obtained from the effective theory, we arrive at the 
identification g — m — 2, which implies that at the bottom of the band the mean occupancy is ^. This may seem 
surprising, as it means that near the bottom of the band there is one particle per two states, a situation that does 
not correspond to neither fermions {g — 1) nor bosons {g = 0) [23|. However, we can argue that this result makes 
indeed sense. Let us first emphasize that we have assumed that the droplet is an ideal gas, which seems to be at odds 
with the fact that the actual particles (i.e. the bosons) are interacting. However, we must remember that under the 
combined action of interactions and rapid rotation, the bosons effectively become hard-core so that their many-body 
wave function has zero interaction energy. This pushes the atoms to orbits where their relative angular momentum 
is equal to 2fi. As a result, there is, on average, one boson per every two angular momentum states: We have N 
particles in the states from ^ to Imax — 2{N — 1), i.e. in 2N — 1 states, which for iV 1 yields two states per 
particle. Therefore, the rapid rotation plus interactions become a "statistical" interaction, which makes the bosons 
behave as if they are non-interacting "super-fermions" obeying Haldane's exclusion statistics ^24. 29, 30] with g — 2. 
Going back to our discussion of the experimental characterization of the Laughlin state, we have thus shown that the 
appearance of the filhng fraction in the spectral weight of Slab{1,^^) can be interpreted as an sign of the peculiar 
statistical aspects of the Laughlin liquid. Therefore, if one measured the spectral weight of Slab{1,uj) and found it 
to be proportional to this would provide a fairly direct evidence for the fact that exclusion statistics is at play in 
the Laughhn state |3(iT31|. 

Another function of experimental interest is the one-body density matrix. To compute it, one needs to find a 
representation for the bosonic field operator at the boundary. To this purpose, it will be convenient to introduce the 
phonon field (piO), related to the density j{6) by de<j>{9) = 2Tij{0). Expanding it in normal modes: 

= ^+Q0+l^l [e^'«6(0 + e~^'%\l)] , (21) 

with [Q, 0o] ~ i. In order to construct an operator with the same properties as the boson field operator, we first 
notice that if a boson is added at the density j{9) becomes j{9) 6{9 — 9o). By direct calculation using (|^ . one 
can check that [j(6'), 0(0')] = ■^5{9 — 9'), which means that j{9) and m4){9) are canonically conjugate to each other. 

Thus the operator that we seek must be proportional to e-*'"'^^*), since e™<^(*°)j(6l)e-^'"'^(^°) = i{9) + 5{9 - 9o). 
In this construction, the boson appears as a soliton or kink in the field (f){9) with topological charge Q =^ +1. 
Fractionally charged excitations can be also constructed, though they are not physical boundary excitations for the 
droplet geometry. They are created by the operator e"*"^^^"-*, which shifts j{9) j{9) + -^5{9 — 9o) and therefore 
describes the creation of a quasi-particle at Zo ~ i?e*^°. One important consistency test, before we proceed any 
further, is to show that the putative boson operator, and its herniitian conjugate, are commuting at different points. 
For 9 9', we have that 

gim0(e)g±im0(e') ^ gTm'[0(e),0(e')]g±irn0(6/')g"n</>(e) 

— (— l)™e**™'^*'^'''e*™'^*'^'' (22) 

and (—1)'" = +1 since m = 2. Indeed, by repeating this calculation with m replaced by i^^^, one can see that 
commuting fields are obtained only when = m is an even integer (if m is odd, the fields anti-commute, and this is 
the situation usually encountered in the FQHE, where the constituents are electrons). Wen has argued |l3 that for 
v there must be more than one phonon branch in the spectrum, otherwise no field operator can be constructed 

and the effective theory is not self-consistent. Below, we shall give explicit examples of how the additional branches 
appear. To close this part of the discussion, it is worth mentioning that the exchange statistics of the excitations 
created by e-"^^"'' is fractional, since e-"^(^) e"*'^^''') = e-''''^sn(e-e')/m ^~i4,{e') ^-i<j,{e) _ rpj^^^ ^^^^^ exchanging two 
boundary quasi-particles we pick up a phase factor equal to ±7r/m = ±7r/2. These excitations are ^-anyons. 
The complete form of the (boundary) field operator in terms of the phonon field is 

= Re'^) = ^\9) = A e-*™['°«+'^(e)]^ (23) 
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where A is a constant that depends on the way the modes at high I > y/mN are excluded from the above sums. 
The value lo = {N — i) is obtained by comparing with the spectral representation of the one-body density matrix, 
Eq. H24|l . Calculation of the one-body density-matrix at the boundary is now possible using (|21|l and (|23|l : 

G{e) ^ (*t(6))^r(0)) e-""'°(^-^') (e-™'^'(»)e™'>(0)). (24) 

This function has two interesting limits (9 3> i/R in both cases). When the temperature is higher than the phonon 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
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e/7i 

FIG. 4: Behavior of the density matrix at the boundary compared with the power-law 1/9^ . The functions have been normalized 
to 1 at = O.Itt. The power-law decay is accurately followed for small angles, whereas for large angles |G(^)| decays more 
slowly since it is periodic on a finite-length boundary. In this respect, notice that |G(^)| has refiection symmetry about 6 = -n. 



level spacing, i.e. if T hv, then 



'ttT 

Gid) = const. X I — 
' Tiv 



sinh 



fT:Te\ 



(25) 



and therefore correlations decay exponentially at large 0, but have power law form at small 9. In the more demanding 
temperature regime where T <C hv, the density matrix reads: 



G{9) ~ const X e 



-imloO 



(26) 



and exhibits an "almost power-law" behavior with 9 (see Fig. 0J| , which is cut off by the finite length of the boundary. 
The last form, Eq. (|26ll . has been numerically shown to be accurate for a system of iV = 36 bosons by Lee and 
Wen (see also Ref. 34 for m = 3 calculations). It is interesting to compare Eq. H26() with the one-body density 
matrix away from the boundary |l35i] . 



G(r,r') = 



-1^ 



z)/2r 



(27) 



Apart from the last term, which is a phase factor, it can be seen that this function decays at large distances as a 
gaussian. This is in contrast to the almost-power-law decay that the same function exhibits near the boundary, where 
|r| « |r'| « R. The differences in behavior are due to the quantum critical fluctuations of the surface modes. 

Nevertheless, it seems difficult that the behavior of density matrix near the boundary, Eq. H26|l . can be experimen- 
tally measured. This is because for a small droplet any measurement of the one-body correlation function would be 
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FIG. 5: Density profile of tfie f — 1 and v = ^ Laugfilin states for A*" = 6 particles. The latter corresponds to boson and the 
former to fermions. The density is normalized to the bulk values po — l/ixP for the u — 1 and po = 1/2-kP for ^' ~ \- Notice 
that the i' = \ state is more extended than the v = 1. The arrows indicate the position of the semiclassical radius R = V mN £. 
Data for the = § state has been taken from Ref. 



dominated by the bulk signal coming from H27|) . However, it turns out that the result (|26|l does have some effect on 
the density profile. To see this, we first notice that Eq. (|26|l can be written using the binomial expansion as 

G{9) = const. X g ^^p^ (28) 
^-^ p\(m — 1) 
p=0 ^ ^ ' 

from which the occupancy of the levels near Imax can be obtained: 

(nil)) oc / deG{9)e'^'^ (29) 
Jo 

with ^ = 0, 1, . . . Thus we find that, for I < Imax, 

(n(0) = const. X ' ^ + ' (30) 

and vanishes for I > Imax, as expected. Therefore, for m = 2, the occupancy near the boundary behaves as (Imax+^—l)- 
But given the limitations of the effective theory, this result will valid for {Imax — I) < VmN. Mitra and MacDonald ^3 
computed {n{l)) numerically for small Laughlin droplets with m = 3, 5, 7 and found good agreement with the general 
expression, Eq. H30|l . within its validity range. Moreover, they also found that (n^l)) exhibits a prominent peak 
followed by smaller oscillations as I decreases from Imax to 0, i.e. as one moves from the boundary to the center of 
the droplet. For large enough droplets, the oscillations eventually damp out as {n{l)) approaches, rather slowly, the 
average value of — . This behavior can be understood on the basis of two facts: i) The existence of a cut-off at zero 
temperature, for / — Imax, above which {n{l)) vanishes according to Eq. (|30|l (e.g. as {Imax + 1 — for m = 2), 
and ii) the fact that the average occupancy equals (= 5 in the present case). These two facts imply that as I 
approaches Imax from below, the occupancy must necessarily decrease below the average |. However, to maintain 
the average occupancy, the particles removed from the neighborhood of Imax must be placed in orbitals with lower 
angular momentum. Since for relatively large droplets, {n{l)) « ^ for I ^ Imax, it seems reasonable to expect that 
{n{l)) exceeds the average by displaying a maximum before it decays to zero as required by Eq. (|30|l . As the orbitals 
with / quanta of angular momentum are located around r; ~ iVl, the peak in {n{l)) translates into a peak in the 
density near the boundary. This can be observed in Fig.|Sl where the density profile of a = i Laughlin droplet with 
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N = 6 particles has been plotted. The density is related to {n{l)) by the following expression, 

Po(r)-XjMO)|^Kr)r (31) 

This result is a consequence of the Laughlin state being an eigenstate of the total angular momentum. When summing 
over I, the occupancy at angular momentum hi is averaged over the neighboring orbitals, and in the present case any 
oscillations displayed by {n{l)) are washed out. The density near the center of the trap reaches the constant value 
of l/(27r^^), which corresponds to the average occupancy i. Nonetheless, a peak in the density near the boundary 
appears and should visible in an experiment where the density profile is measured. In this respect, it was shown in 
Rcf. 21] using Eq. (|31|l that by turning the trapping potential off and allowing the atoms to expand freely, the density 
profile of a vortex liquid expands self-similarly. This fact can be used to perform very accurate measurements of the 
density profile, which should reveal a characteristic peak near the boundary for the Laughlin state (and possibly for 
other vortex liquids), and which is not expected for Bose-condensed systems. 

For comparison purposes, in Fig. [S] we have also plotted the density profile for the fermionic ly ^ 1 Laughlin state. 
This is a state where particles are non-interacting, and the occupancy is known exactly: {n(l)) — 1 for / < Imax = N—1, 
and zero for / > Imax- This is again in agreement with (|30l) . which predicts n{l) = const., for I < Imax (this constant 
cannot be fixed by the effective theory). In this case {n{l)) vanishes abruptly for I = Imax + 1- Therefore, we do not 
expect a peak in the density near the boundary. We conclude that the peak is a consequence of the strong correlations 
built into the Laughlin state. It is indeed a generic property of finite droplets of the Laughlin liquid, and in a more 
general framework is a particular example of the generalized Luttinger's theorem p5| discussed by Haldane in Ref . 

We end this section by emphasizing the fundamental difference between the boson occupancy, {n(l)), and the 
occupancy for the ideal "super-fermion" gas, N{1), introduced in our discussion of the density response of the Laughhn 
droplet. The latter is a convenient tool, which can be also used to discuss some other properties of the Laughlin 
liquid, e.g. the thermodynamics of the boundary excitations or the Hall conductance p3l| . The boson occupancy 
(n(Z)) contains in addition correlation effects between the actual particles, which lead, for instance, to the almost 
power-law behavior of G{9). A mathematically rigorous discussion of these issues can be found in Ref. 13 II. 



IV. OTHER STABLE STATES OF SCALAR BOSONS 



The Laughlin state described in the previous section has the property of being an exact eigenstate of the Hamilto- 
nian: 

N N 

H = Y,Hj+g2d '^(r-r,). (32) 

i=l i<j=l 

More precisely, ijHJ) is the state with the lowest energy for L = fimN{N — l)/2. On the other hand, the states that we 
will consider in this section are not exact eigenstates of H32|l . They have been obtained numerically using two different 
types of exact diagonalization set-ups. If one is only interested in the bulk properties of vortex liquids, and taking 
into account that exact diagonalization is only feasible with relatively small particle numbers, the Hamiltonian can be 
efficiently diagonalized on a compact manifold, like a torus or a sphere @ . Cooper et al. used the toroidal geometry 
to study the stability of the Abrikosov lattice as the filling fraction v — N/Ny {Ny being the number of vortices) is 
varied. They found that the lattice melts for ^ ~ 6 and that for smaller filling fractions a rotating bosonic system 
exhibits a series of gapped states {homogeneous vortex liquids) when ly belongs to the sequence i, 1, |, 2, . . . , | . . . The 
gap is a signature of the incompressiblity of the vortex liquid. The = ^ state was found to correspond to the Laughlin 
state discussed above, while the remaining states were shown to have good overlap with the Moore-Read state (for 
= 1, see below) and other wave functions of the class introduced by Read and Rezayi 27j. At the same time, the 
Laughlin and Moore-Read states had been previously found in another type of exact diagonalization studies [l(l llll | 
where the bosons are confined by a harmonic potential. This is a situation that is closer to the one experimentally 
relevant, and in which the vortex liquids are inhomogeneous. In this set-up. Cooper and Wilkin [loi | found a series of 
stable states at magic values of the total angular momentum, L. The series terminates at the f = ^-Laughlin state, 
and many of remaining the stable states have good overlap with wave functions constructed from the following ansatz: 

<i>'^\zi,...,ZN)^VLLL [$i(zi,...,Z7v)$^"Hri,...,r7v)] , (33) 

where $1 = ^^^^ denotes the Jastrow factor, ni<j(-^i ~ -^j)' ^ Slater determinant representing the state 

of N fictitious fermions, the so-called "composite fermions" (CF's), distributed over n Landau levels. The operator 
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(b) 
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FIG. 6: A compact state with = 9 composite fermions (a and b), which corresponds to the integer sequence [4, 3, 2]. In (b) 
we have represented the effect of the confinement, which is to lift the degeneracy of the otherwise degenerate levels within every 
CF Landau level. Surface excitations can be pictured as low-energy particle-hole excitations (b) within each Landau level. For 
these excitations the composite fermions can be weakly interacting, even if the interactions could be neglected in the ground 
state. 



Vlli, is needed to project the ansatz onto the LLL (for details of how this projection should be carried see Ref.0 and 
Ref . l30|) . Furthermore, Cooper and Wilkin found that the stable states are "compact", a term coined by Jain and 
Kawamura to refer to those states where Ni CF's occupy the lowest angular momentum orbitals in the i-th CF 
Landau level (« = 0, 1, . . . , n — 1, so that X^T^o — Fig. EJ- Therefore, a compact state is characterized by 

the set of integers: [A^o, Ni, N2, . . . , A^n- il, a nd it has an energy equal to Tiujcf J2^=q + where ujqf — ^^cf{N) 
is the effective CF Landau level spacing 1l(ll28l |. However, it is necessary to stress that not all compact CF states are 
stable 28]. In the large A^ limit, the CF ansatz yields homogeneous vortex liquids whose filling fraction falls into the 
bosonic Jain sequence i.e. = n/{n -1- 1) or = (n + l)/n, for an integer n > 1. Except for the v = i-Laughlin 
state, this sequence does not match the one found by Cooper et al. using the toroidal geometry, and this has led 
to some confusion. The situation has been recently clarified by Regnault and Jolicoeur _6J who performed exact 
diagonalizations on another compact geometry, namely a sphere, and checked for the convergence of the spectral gap 
with the system size. They found the existence of stable homogeneous vortex liquids for = i, |, |, |, |, falhng into 
to the Jain sequence, and ior v = 1, which corresponds to the Moore-Read state. Therefore, below we shall consider 
the surface modes of some of the states in the Jain sequence as well as the Moore-Read state. 

To study the surface waves of the vortex li quid s whose wave functions can be obtained from the CF ansatz, Eq. H33|l . 
we will employ the parton construction [T9I l37j. Usually, this method assumes that the number of particles, A^, is 
macroscopically large so that the inhomogeneity of the system can be neglected. For these states, one can safely 
assume that there are n CF Landau levels filled by N/n composite fermions. As pointed out above, the resulting 
state is characterized by a filling fraction v = n/{n + 1). However, when dealing with a mesoscopic droplet, this is not 
necessarily true and the compact states exhibiting good overlap with the stable states usually have unequal number 
of CF in different Landau levels. Generalizing the parton construction for this situation is not difficult, the only 
limitation being imposed by the subsequent use of the effective low-energy theory introduced in Sect. Illll to describe 
of the surface waves of a given CF Landau level, i. This requires that the CF level occupancy Ni ^ 1, which is 
not always fulfilled. Nevertheless, it was found [Sg that having some of the higher CF Landau levels occupied with 
Ni ^ 1 CF's corresponds to quasi-particle excitations over a state of the Jain hierarchy. Therefore, these states can 
be characterized by an effective number of CF levels to which the effective theory can be applied. In the discussion 
that follows, we denote this number by p (< n). In other words, this integer, which yields the effective filling fraction 
f = p/{p + 1), can be associated with the number of CF Landau levels with occupancy A^i ^ 1 . 

Keeping in mind the above caveats, we can regard the CF ansatz H33|) as describing a bound state, the fundamental 
boson, made up of two kinds of fermions (called "partons") such that the many-body wave function (we drop T'lll 
for notational simplicity) reads: 

<\zi,...,ZA.) = ci>i(z(^\...,z«)$(")(zf\...,z(,^)) . (34) 

-l z. —Z . —Zi > 
I. • ' J i=l,...,N 

Thus a boson is a type-1 parton in the (fermion) LLL, which carries a "charge" qi = p/{p + 1), bound to a second 
parton (type-2), which carries charge (72 — l/(p+l) so that qi+q2 = 1. Type-1 has one branch of boundary excitations, 
whereas type-2 partons contribute with p < n branches corresponding to the p CF levels with A'i ^ 1 (to see how this 
comes about, repeat the steps that lead to the effective theory, now setting m = 1 for type-1 partons, and considering 
p different branches with m = 1 for type-2 partons). However, not all of these p + 1 branches are independent. In 
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fact, the bound state condition implies that the density of type-1 partons cannot fluctuate independently of that of 
type-2 partons. Therefore, we must demand that all physical operators do not excite density fluctuations where both 
partons behave independently. Mathematically, if O is a physical operator, then [pc, O] — 0, with = pQ — J2^=i Pi: 
where po describes the surface modes of type-1 partons, and pi for i ^ 1, . . . ,p those of type-2. These operators obey 
the current algebra: 

[p,{l),p,il')] = l'5i+i,.oS,,,. (35) 

As mentioned previously, there are p independent branches, which are obtained by requiring the currents to be physical 
operators, in the sense introduced above. Thus it is found that one of the independent currents corresponds to the 
total current jo = qipo + 92 Y^^=i Pi, the remaining p — 1 are given by ji — '^ijPj (* — 1: ■ • -P ~ 1); where 

Uij are p — 1 orthogonal vectors such that X]j=i ~ 0. 

To illustrate how the parton construction works in practice, we first notice that one can readily recover the results 
of Sect, mil for the surface waves of the v — ^ Laughlin state, which in the present context corresponds to setting 
p = n = 1. However, a less straightforward exercise is provided by a state where p — 2 (and v ~ ^ effectively). Using 
the parton construction, we find that there are two independent branches of surface modes, described by the currents 
jo = (2po + Pi + P2) /3, and ji ~ pi — p2, respectively, which obey 

Ml), Ml')] = ll'Si+,^0, (36) 

Ml), Ml')] = 2l'Si+i,,o- (37) 

As a check, it is worth trying to obtain these results from a different point of view. Let us first recall that the v ~ ^ 
state can be obtained by means of the hierarchical construction 38] where the state at v — 2/i can be thought of as 
two component vortex liquid, containing one liquid with 1^0 = ^, which corresponds to the parent Laughlin state, and 
another one with 1/1 = ^, which corresponds to a Laughlin state of quasi-particles over the parent state of bosons. 
Notice that the sum of filling fractions + vi — v. If, for the moment, we neglect any interactions between the two 
components of the liquid, the effective theory can be used to describe each one separately. This leads to the following 



current algebra: 

[ao(0,'To(/')] = ^^'W,o, (38) 

[ai(0,ai(?')] = ^/''5z+r,o, (39) 
6 

where ctq describes the boundary phonons of the = ^ liquid and o\ those of the v\ — ^ liquid. 



To find the relationship between the parton currents and the currents ctq smd a\ , we resort to the interpretation of 
jo as the total current. This implies that jo = gq^g\. The other current must be an independent linear combination, 
which is readily found to be ji — — 3cri. This, quite direct, identification of the currents helps to confirm the 
results obtained from the parton construction. Thus we conclude that, just as we briefly remarked in Sect. IIIII for 
states other than the Laughlin states there are several branches of surface modes. In particular, for the compact CF 
states considered above, there can be p > 1 branches. 

Let us finally turn our attention to the dynamic structure factor for the states considered above. As it was pointed 
out previously, these states are not exact ground states of the Hamiltonian, which means that we cannot use the wave 
function approach of Sect. ITTll to get much insight into the excitation energy of the surface modes. Furthermore, the 
non-interacting CF picture may only hold to describe the ground states but not their low-lying excitations. In other 
words, some small residual interactions are always expected. However, if we start by assuming that the non-interacting 
CF holds, then the energy of the surface modes will be just the confinement energy, and therefore the Hamiltonian 
reads: 

fl"o = ?i7ruV / dejl{9). (40) 

In the above expression we have normalized the p independent currents so that [Jq,(/), Jg(r)] = I' 5a.i35i+i' ,q (notice 
that this normalization differs from the one used above for jo and ji, i = 1, . . . , n — 1). In the absence of interactions, 
the p phonon branches are all degenerate. However, since the potential 5vlab{Q) (see appendix^ couples only to 
jo(0) — y/T'Jo{0) (where v =p/{p+ 1)), Slab{1,^) will have the form (|17|l . displaying a single peak at a; = uj±l, but 
with pI [p -f 1) replacing 1/m or, in other words, with total spectral weight proportional to p/ {p + 1). 
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The spectrum is modified, however, if the residual interactions could not be neglected. In this case, one has to 
consider the effect of adding to the Hamiltonian terms like: 

H'^n^Y, / ^Vo,pJo.{0)Jp{e), (41) 

where Au^/j is a real symmetric matrix. One could also have included terms containing derivatives of the currents, such 
like gap J d9 Ja{6)deJi3{0). Other terms containing more than two current operators are also possible, and describe 
inelastic scattering processes among the surface phonons. Being irrelevant in the renormalization-group sense, these 
terms yield corrections to the energy of the surface modes of the order of l{l/\/mNy , where the index i equals the 
order of the current derivatives plus the number of current operators minus two (e.g. i = 1 in the previous example 
containing two currents and one derivative). These corrections can be important for modes with I ^ \/N . To the 
lowest order in this case, however, the "normal" (surface) modes can be obtained by diagonalizing the matrix 

TiaP = V Saf} + ^Vap. (42) 

Nevertheless, the overall shape of the spectrum is constrained by the generalized Kohn's theorem discussed in ap- 
pendixEl This theorem is a consequence of the decoupling of the center of mass motion from other degrees of freedom 
in harmonically trapped systems. Since we have assumed that this is the case throughout, it implies that the energy 
(in the laboratory frame) of the / = 1 mode is given by the hare trap frequency lo±^ . Let us begin by considering the case 
where there is a single phonon branch, like in the the Moore-Read state to be discussed below. Kohn's theorem implies 
that Lo{l = 1) = uj±. Any corrections arising from interactions will necessarily have to vanish for / = 1. For instance, 
we can think of a general form for the dispersion (in the laboratory frame) like = LOi_l + ^X^S" ^'^inti^ ~ 1)% 
where Suj^l^^ ^ lu±. In the case where there is more than one branch, the situation is slightly more complicated. 
Then, in principle, after diagonalizing the matrix Tiap, not all the eigenvalues will be degenerate, which implies that 
for I — I, to this order, there can be several modes with distinct energy. However, the energy of, at least, one of the 
modes (the Kohn mode, associated with the center of mass motion) is still fixed to uj± (in the laboratory frame). 
Furthermore, an analysis of the oscillator strength based on Kohn's theorem (see appendix IXI and below) shows that 
the Kohn mode exhausts the all the oscillator strength available, which in turn forces all the modes to be degenerate 
at / = 1. Thus, the general form of the dispersion of the normal modes is a generalization of the single-branch case: 
uj{l, a) = uj±l + lJ2t^ ^^inti'^)i^ ~ 1)% where a = 1, . . . ,p is the branch index. This form assumes the frequencies to 
be real and thus neglects any imaginary parts arising from inelastic coupling of surface modes. However, this fact can 
be readily accounted for, although we expect this linewidth to be small for the low- lying modes, and zero for 1 = 1. 

The implications of the previous discussion for the dynamic structure factor can be easily extracted. After expressing 
the total current jo — y^Jo in terms of the normal modes, we find that Slab{1,(^) should exhibit the following 
structure 

lwaS{uj-u{l,a)) 

Q — 1 

From the current algebra for jo(Oi follows the following sum rule: 

j2wa = ^. (44) 

According to the above considerations, for Z = 1 all the branches must be degenerate: w(l, a) = uj±, and hence there 
is single peak of spectral weight proportional to p/{p + 1). For I > 1 more than one (but at most p) peak can exist, 
and given that for repulsive interactions the energy of the Kohn mode is a lower bound of the surface phonon energy, 
all the peaks should appear at frequencies higher than (or equal to) uj±l. However, since the interactions are weak, 
it may well happen that the shifts and splittings are very small and hard to resolve experimentally, and this would 
produce a single peak again. The total spectral weight associated with the surface modes must be again proportional 
to p/{p + 1), for integer p > I. However, if the different peaks for ^ > 1 could be resolved, the individual spectral 
weight of each peak would be determined by the detailed form of the interactions between modes, which is not fixed 
by the present effective field theory. 

We close this section by considering the surface modes of the Pfaffian or Moore-Read (MR) state. This is a state 
that cannot be constructed from the the composite- fcrmion ansatz, Eq. H33|l . As it was mentioned at the beginning of 
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the section, for v = 1 exact diagonalization studies on the torus [Bl and the sphere have found an incompressible 
state exhibiting a large overlap with the MR state, whose wave function reads: 



$^^(zi,...,zjv)=^ 



21 - Z2 ZN-1 - ZN 



N 



n ~ 

i<j=l 



(45) 



where A[. . .] means that the bracketed product must be anti-symmetrized as follows (N is assumed to be even): 



A 



1 



Zl - Z2 ZN~1 - ZN 



1 



2JV/2(7v/2)! 



1 



1 



(46) 



In the previous expression P is a permutation of 1,2, . . . ,iV and (—1)^ denotes its signature. For instance, if we 
consider iV = 4 particles. 



A 



1 



1 



Zl - Z2 Z3- Z4 



1 



1 



1 



1 



Zl - Z2 Z3 - Z4 Zl - Z3 Z2 - Zi Zl - Z4 Z2- Z3 



(47) 



Alternatively, the MR state can be written in the following way |^ ^| : 



MR 



(zi, . . .,zn) = S 



{i<])eA 



ik<l)eB 



(48) 



where S symmetrizes the product of the two Laughlin i states over all possible partitions of N particles into two 
disjoint subsets A and B having N/2 particles each. The second form, Eq. (|45|l . shows perhaps more explicity that 
the lower angular momentum of the MR is attained by having N{N — l)/2 — N{N/2 — l)/2 = {N/2)'^ pairs of bosons 
in orbits with zero relative angular momentum. This is to be contrasted with the Laughlin state, Eq. (jS)), where all 
boson pairs have relative angular momentum equal to 2h. This means that in the MR some particles must necessarily 
interact, in order to reduce the total angular momentum, which for this state equals N{N — 2)/2 for N even and 
L = (TV — 1)^/2 for N odd. Wilkin and Gunn found in Ref. [Tl! that this state is also a good candidate to describe the 
stable state found for rotating bosons in a harmonic trap at L = N{N — 2)/2 (for even N). 

Surprisingly, the MR state becomes an exact eigenstate with zero interaction energy for a model where the particles 
interact by means of a repulsive three-body potential (53 > 0) |23| 



V3=g3 E ^''Hr.-rfe)(5(2)(r,-r,). 

i<j<k 



(49) 



This follows from the property that MR vanishes when any three-particles come together, and has implications for the 
stability of the state with respect to three-body losses. We have previously found that the Laughlin state was a zero 
energy eigenstate of a two-body interaction potential. This is because its wave function vanishes as any two-particles 
approach each other, hence the probability for three particles to be at the same point must be necessarily zero as 
well. Therefore, for the Laughlin and the MR states the three-body recombination rate 43] should be zero, rendering 
them very stable. However, the MR state is only a good approximation to the true ground state of a boson system 
interacting with two-body potentials, which implies that in the actual system the recombination rate should be small 
but not zero. 

However, the existence of a Hamiltonian for which the MR state has zero energy allows to perform a microscopic 
analysis of the surface excitations formally analogous to the one carried out for the Laughlin state. Thus it is found 
that the symmetric polynomials si — X]i3=i ^1 describe, also in this case, low-lying excited states of the Hamiltonian 
V3 introduced above: The states obtained multiplying the MR state by s; have zero interaction energy, and therefore 
are degenerate with the MR state, but adding the confinement term, vL, this degeneracy is lifted. We can identify 
the surface modes of the MR state with the polynomials s/, exactly as it was done for the Laughlin state. But it has 
been shown by several authors "2^, 'S] that these excitations do not exhaust the spectrum of boundary modes of the 
MR state. Besides the chiral phonons, which have bosonic character, fermionic excitations that do not carry charge 
quantum numbers (i.e. they are described by a Majorana fermion) also exist. Lack of charge quantum numbers means 
that it is not expected that these excitations will be directly created by weakly deforming the confining potential (it 
can happen that the fermionic neutral and the bosonic charged excitations are somehow coupled, but the form of the 
interaction seems not to be easy to write H^). Thus we will focus only on the single branch of "charged" surface 
phonons. By the same arguments employed with other states of scalar bosons, the phonon energy is given (in the 
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laboratory frame) by w = uj±l. But since the MR state is indeed an approximation, we cannot exclude the possibility 
that the phonon dispersion will receive corrections for I > 1 (for I = 1 the energy is fixed by Kohn's theorem, as 
discussed above). However, this does not affect our conclusion that, for the MR state, the dynamic structure factor 
should display a single surface phonon peak, whose energy is approximately linear with the angular momentum of the 
mode. One can repeat the steps that lead to the effective field theory for a droplet of MR liquid, and obtain that the 
spectral weight of the surface phonon in Slab{1,'-^) is again proportional to l/m with m = 1. 

V. VECTOR BOSONS 

Let us now take up bosons with internal (i.e. hypcrfine) degrees of freedom. We call them vector bosons because 
we will describe them by a field operator that is a tuple of n-fields, 

*t(^)^(.l.t(^),...,vl;t^_^(^)), (50) 

and which transforms as a vector under SU (n) transformations. This does not necessarily mean that the Hamiltonian 
has this symmetry, but we will find that under certain conditions the effective theory describing the surface modes 
does. One can also regard the internal degree of freedom as a spin or pseudo-spin, depending on the context. Thus 
an ri-component vector boson corresponds to a boson of spin S ^ (n — l)/2. Again, this does not mean that the 
Hamiltonian is in general rotationally invariant (see below). Scalar bosons (i.e. S = 0) are obtained by letting n = 1. 
For n = 2 we obtain spin-1/2 bosons, which can be used to describe a mixture of two hyperfine states of the same 
isotope (as in the experiment reported in Ref. |^). If the atoms are trapped by purely optical means, one can have 
in the same trap |5Clj| the three magnetic sublevels (i.e. mp = —1,0, +1) of an isotope of total spin F — 1 IsBI, and 
this case will be treated by setting n — One also can imagine experiments with atoms of higher S (F). In general, 
however, the interaction will depend on the hyperfine state of the colliding atoms. Thus for S = 1/2, it can take the 
generic form 

V^, = [go + giiA.S'tSJ + AySfS] + A^S^S^)] 5(r, - r,), (51) 

for any two particles i and j. For spin-1 bosons in an optical trap rotational invariance implies that |49| | 

- [go+9lS^■S,]S{n-r,), (52) 

Higher spin will involve higher powers of (Si • S^) j49j |. Typically one has l^ij ^ go, while the sign of gi, which 
determines the ferromagnetic (i.e. gi < 0, like in ^^Rb) or anti-ferromagnetic (i.e. gi > 0, like ^^Na) character of the 
spin-interaction, depends on the atom species. 

We have seen in previous sections that as the amount of angular momentum deposited in the system is increased, 
all atoms tend to avoid each other. When the atoms have internal degrees of freedom, a cloud containing an equal 
number of atoms in each internal state will go to a singlet state with zero interaction energy at sufficiently high 
angular momentum. This possibility has been recently considered by a number of authors for spin-1 [2lll2^ as well as 
for arbitrary spin bosons 18]. The parton construction can be readily applied to study these states. We will consider 
that the interaction is dominated by the term proportional to go, but including a weak (as it is the usually the case) 
spin-dependent scattering term makes no difference. We stress that this is because, as long as the interaction has 
zero-range, these states have zero interaction energy because particles completely avoid each other independently of 
their internal state. 

To see how the parton construction allows us to describe singlet states, we begin by writing the components of field 
operator, Eq. H50I) . as the product of two field operators of the fermionic partons, i.e. 

^iiz)^4iz)4jz), (53) 

for a — 0, . . . ,n — 1. In terms of wave functions this implies that the boson wave function reads 

$B(ziai, . . . , znUn) = '^i{zl, . . . , zlf)^^"\zfai, z^un) 

where $i = ni<j(-^i^ ~ ^j) denotes the Jastrow factor, which describes the state of type-1 partons, whereas <I>("^ is 
a Slater determinant of N type-2 partons in internal states ai, . . . , un- In this context we are interested in states 
where all type-2 partons are in the LLL. The situation differs from the case of scalar bosons because now type-2 
partons carry the internal degrees of freedom, and therefore states other than the Laughlin state will appear for which 
the projection onto the LLL is unnecessary. It is interesting to point out that it is the Jastrow factor that ensures 



(54) 
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that the resulting wave function has zero interaction energy, while the kinetic energy is quenched to Nhuj± plus 
the confinement energy. In other words, the factor $1 makes the wave function vanish whenever any two-particles 
approach each other, independently of their internal quantum state. For the same reason, the probability for three 
particles to approach each other is zero, which in turn implies that the three-body recombination rate vanishes: These 
states should be particularly long-lived. 

To make things more concrete, we specialize our discussion to n = 2 components, i.e. S — ^ bosons. We consider 
this case first because it is simpler but already illustrates the important points, and also because it is relevant to the 
situation where only two hyperfine states are allowed in the same trap j48l |. We will discuss the generalization to 
arbitrary n at the end of this section. Setting n = 2, we assume that a — corresponds to spin up (t) and a — 1 
to spin down (|). As mentioned previously, we shall consider only singlet ground states. Therefore, the number of 
particles, TV, must be divisible by n = 2, i.e. iV| = iVj^ = N/2 . Nevertheless, the latter condition does not suffice for 
to be a singlet but only ensures that Sz\^b) = 0. In addition it is needed that 



5-|$r) =0, 



(55) 



where S = X]j=i = '}lj=i{^x ~ ^^y) is the operator that decreases the total spin. As the spin index is carried 
by type-2 partons, the above condition must be indeed met by the parton wave function <I>(^-'. Upon filling the LLL 
with iV/2 fermions of both spin species, the parton wave function takes the form 



N/2-1 

1=0 



(56) 



where |0) is the zero particle state. It is not hard to see that this state obeys the condition H55|) since the operator 
= X]/!=^ V'2x (0V'2t (0 annihilates it. However, in order to obtain the boson wave function, we need the spatial 
dependence of H56I) . This is just the product of two Slater determinants of N/2 fermions in the LLL, one for each spin 
orientation, which have the form of the Jastrow wave functions that we have been using throughout. Hence, using 
Eq. one finds 



$221 (^1, ■ • . ,Zw/2; Wi, . . ■,Wn/2) 



N/2 N/2 

n ~ ^jfi^^ - n 

i<j=l ij=l 



(57) 



where and Wi (z = 1, . . . , A'^/2) denote the positions of the spin up and down bosons, respectively. Notice that this 
wave function seems to have a rather peculiar dependence on Zi and Wi, which is not fully symmetric in these variables, 
in spite of being a bosonic wave function. However, one must be careful since we are dealing with spinful bosons. 
Indeed, the above form is what the singlet condition requires. As the reader can easily check for two bosons making a 
singlet, the anti-symmetry of the spin part of the wave function requires an anti-symmetric spatial wave function. The 
appropriate generalization of this observation to many particles is given by (|57(l . The <I'22i wave function is indeed 
one of the class introduced by Halperin 0,113 to describe non- fully spin-polarized quantum Hall states. In Ref. I23I 
this state was characterized by the affine algebra sm(3)i, i.e the symmetry of the two-dimensional (conformal) field 
theory for which the wave function can be obtained as a correlation function (see the appendix for a simple example, 
and below). 

After explaining how to construct some of the vector boson states from partons, we proceed with the description 
of their surface modes. Indeed, all we need to do is to re-interpret what was done in the previous section for the 
n = p = 2 state. Now jo must be interpreted as the charge current, which we denote as jc in what follows, whereas 
ji/2 measures the density of the third component of the spin at the boundary, and we denote it as j^. Therefore, 



MO) 



Po{0) 



(PT +Pi)(^) 



(58) 
(59) 



A difference with respect to the case of scalar bosons is that, in this case, jc and are decoupled, as long as the 
trapping potential is spin-independent. The symmetry of the problem is C/(l) x SU{2), where where U (1) is associated 
with the total particle number (i.e. the charge), and therefore with jc, whereas SU{2) is related to the internal degree 



of freedom, and therefore to the spin current, j'^ 
which also includes the currents 



Indeed, the spin current is part of the generators of a larger algebra, 



V'2|'02i -t- V'2J.'^2T 
^2TV-'2i - V'2iV'2T 



(60) 
(61) 



18 



These three currents can be expressed in a more compact way as jg{9) = 5 X]q /3 '^'ai^)^a0^pi.^)^ ^^"^ they obey the 
foUowing algebra (a, b,c = 1, 2, 3): 

[Jsil),j'in] = |'5''V;',o + ^efm + 1% (62) 

known as sm(2)i Kac-Moody algebra. 

If decoupled, the charge, jc, and spin, , currents describe waves that propagate with the same frequency v. This 
can be shown using the wave function approach introduced in Sect. lIIll In this case, each spin orientation is related to 
different set of symmetric polynomials, U = Y^fl^ z\, ui — J^iHi "^i^h the same excitation energy, hvl. The charge 
modes corresponding to the combination U + ui and the spin modes to ti — ui, are degenerate. Thus we conclude 
that the dynamics of these modes is dictated by the confining potential (which we have assumed spin-independent). 
Just like for the Laughlin state, this is a consequence of the wave function $221 being an exact eigenstate of the 
Hamiltonian with zero interaction energy. As far as the observable consequences are concerned, it is now possible 
to consider two dynamic structure factors, given by the correlation functions {jc{0,t)jc{0)) and {j^{6,t)j^{0)), and 
related to deformation potentials that couple to the total density, jc, or the spin density, j^. In both cases, the 
observable spectrum will exhibit a single peak at w = uj±l. The difference will be in the spectral weight of the peaks, 
which for the charge peak will be proportional to = 2/3, whereas the spin peak to 1/2. We finally point out that 
if the confining potential is spin-dependent, then the separation into a charge and a spin density current is no longer 
convenient. In other words, the currents jc and do not correspond to the normal modes of the system any more. 
In such a case, for instance, the spectrum for a trapping potential deformation that couples to the total charge will 
display two peaks at two different frequencies to — ujijj_l and uj = uJi_\_l (we assume that W|.x_l — ^ ^ The total 
spectral weight of the two peaks, however, will be again proportional to 2/3, each peak contributing 1/3. 

As a final demonstration of the internal consistency of these constructions, we will rederive the wave function $221 
using the formalism explained in the appendix, and which relates the wave function to the correlation functions of 
the effective theory for surface modes. The key to this approach is to identify the boson operator in the LLL with the 
vertex operator (see appendix for a definition) that describes the boson operator at the boundary. To this purpose, 
let us introduce the phonon fields for the partons, in analogy to how it was defined in Sect. Hill i.e. 



MO) = 27: Pole), (63) 
deMO) = 2^p„(0), (64) 

with a =1, 1. The the parton construction dictates that the boson field operator at the boundary takes the form [T9l | 

^t^(0) =e-*('^"+'^°)(^), (65) 

which is nothing but Eq. I|53(l with the parton fields written as exponentials of their corresponding phonon fields (cf. 
Sect. IIIl)l . However, the phonon fields (po, (f>^ , (pi, just as the currents po,p^,pi, are not independent. Only jc and 
are independent. But the following expressions: 

Po + PT = Ijc + j!, (66) 
Po + Pi = ^jc~js, (67) 

help us to relate the phonon field combinations that enter in ^'J^, to the charge and spin phonon fields, defined by the 
following equations: 

de^cie) - 2n^jc{e), (68) 

dgipsie) = 2TrV2jU0). (69) 
These fields have been so normalized that their vertex operators have the two-point correlation functions: 

{Vli{z)Vl,{0)) = (: e"^'(^) : : e"''^^'*" :) = ^, (70) 



where i,j = c, s. The field operators read 



^^^^(z) = :e~<^^^+7f^0('^) : , (71) 
^l^iz) = :e-KVi'^=-75'^0(^") : . (72) 
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We are now ready to compute the wave function. It takes the form 
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(73) 



This is just the anti-analytic version of the Halperin 221 state, Eq. H57|) . In the above expressions R is the radius of 
the disk of positive charge background with density po = vjnP' = 2/(37r£^) and total charge equal to N . It is also 
interesting to notice that, in the above construction the wave function splits into two factors, one involving only charge 
(}Pc) and other involving only spin (}Ps) operators. This decomposition can be understood as due to the separation of 
the fundamental boson into two kinds of excitations: a spinon, which carries the spin, and a holon, which carries the 
"charge" . 

Finally, the generalization to arbitrary n is not hard to work out. For a number of particles divisible by n, the 
singlet SU (n) wave function reads: 



K J a—0 



,i-l N/n 

n n - 4f 

a=Q i<j=l 



N/n 

n n « 



(74) 



where a, /3 = 0, . . . , ti — 1 (This state is the su{n + l)i state of Ref. [23j). Just as we did above, one can construct 
the surface modes from the symmetric polynomials, sj"-* — 'Ylif/i (•^1"^) ' ^^ich will have excitation energy (in the 
rotating frame) hval, with vq = vi = . . . = Vn-i = v = u!± — Q if the trapping potential is spin- independent. The 
parton construction yields n independent currents, one of them, jc = {npo + J27=i + 1)' corresponding to the 

total current and is therefore associated with the charge U{1) symmetry. The remaining n — 1 currents are part of a 
su{n)i algebra that describes the dynamics of the internal degree of freedom. The dynamic structure factor related to 
the {jc{0, t)jc{0)) correlator exhibits a single peak dX oj — uj±l, whose intensity is proportional to i' — n/{n + l). If the 
trapping potential is spin-dependent, however, this peak splits into as many peaks as distinct values of the frequency 
Wct_L exist, but the total intensity remains proportional to n/{n + 1) each spin component contributing l/(ri + 1). 



VI. CONCLUSIONS AND OUTLOOK 



In this paper we have described the surface modes of ultra-cold atomic clouds with a very large number of vortices. 
These systems are known as vortex liquids, and are expected to be observed in the critical rotation regime where 
the rotation frequency approaches the trap frequency. The surface waves of vortex liquids present several important 
differences with respect to the surface modes of a Bose condensate. For an axially s ymm etric trap of the type 
considered above, the frequency dispersion of the surface modes of a BEC is found to be |49| : 

LU^LU^y/lTl (75) 

where I — 0, ±1, ±2, . . . On the other hand, for a spin-independent trap, a vortex liquid exhibits surface modes whose 
dispersion is well approximated by 

u) = u;_lI, (76) 

where / = 1,2,3,... i.e. a positive integer, which measures by how many quanta is increased the total angular 
momentum of the vortex liquid. We also found other (approximate) excitations, which involve promoting particles to 
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higher Landau levels, and which have the same energy as the surface modes in the laboratory frame. However, these 
excitations decrease the total angular momentum by Z = 1, 2, . . . quanta. However, we must stress that the surface 
modes are the true low-lying excitations of the vortex liquids as they dominate their low-temperature properties (for 
T <^ hQ) since it is their energy in the rotating frame, not that in the laboratory frame, that enters the Boltzmann 
factor. Furthermore, we argue that for T ^ h{LU± — Q), their intensity in the dynamic structure factor should be 
much larger. 

The chirality of the surface modes, namely that I is restricted to positive integers, is a consequence of the dramatic 
effect that rapid rotation has on a cloud of ultra-cold atoms: Besides enhancing quantum fluctuations, which are 
finally responsible for the destruction of the condensate and the melting of the Abrikosov lattice, fast rotation, when 
regarded from the rotating reference frame, has the same effect as a strong magnetic field on a system of charged 
particles moving in two dimensions. In both situations, time-reversal symmetry is broken, which leads to the chirality 
of surface modes. We find, by comparing with the angular momentum and excitation energy of bulk excitations of the 
Laughlin liquid, that the surface waves are the low-lying modes of the system. Indeed, quasi-holes and quasi-particles 
are bulk excitations that change the total angular momentum by ^ hN, whereas the surface modes do it just by 
I ^ V mN quanta. We expect that this property also holds for the more complicated vortex liquids described in 
Sect. IIVI For these states, we find that there exist more than one branch of surface phonons. This result is obtained 
by relying on the composite-fermion ansatz used in Ref . Hol and first introduced for fermions in Ref.l28i In this picture, 
a bound state of a boson with a vortex is an object of fermionic statistics, known as composite fermion, which occupies 
n "renormalized" Landau levels. The parton construction 19] used in our calculations in a concrete implementation 
to study surface excitations of the composite-fermion idea. 

We propose to probe the surface excitations by adding a weak time-dependent part to the trapping potential 
that imparts small amounts of angular momentum to the system. This should excite surface waves in droplets of 
vortex liquid. We have quantified the response of the system in terms of the dynamic structure factor Slab{1,^)- 
This function is related to the density correlation-function at the boundary. If the residual interactions between the 
composite fermions can be neglected, S'las(^w) should exhibit a single peak at the frequency given by Eq. ifTHll . 
However, in the case where residual interactions are not negligible, several peaks may appear for I > 1. We also find 
that the spectral weight of the surface modes in Slab{1, w), is proportional to p/{p -\- 1), where p is a positive integer. 
This means that if the system is gradually driven to states of higher angular momentum, going through different 
vortex liquids at every stage, the integer p characterizing the peak intensity of the surface modes, should decrease 
down to p = 1, which corresponds to the system being in Laughlin state with v = \ (higher angular momentum states 
can be imagined but this would be experimentally much harder). In other words, the total intensity in Slab{It^) of 
the surface modes at a given I should decrease as the Laughlin i state is approached, at constant temperature. 

Finally, we have also found chiral surface waves in vortex liquids of vector (i.e. spinful) bosons. In this case, the 
parton construction, used to study the surface waves of scalar bosons, can be also successfully applied. If the trapping 
potential is spin independent, the experimental signatures are similar to those for states of scalar bosons with p equal 
to the number of components of the boson field (2 for spin i bosons, 3 for S = I, 25* -I- 1 spin-S bosons). However, 
in this case the surface-phonon peak is always to be found for lo = where Z > 0. In spin-dependent traps, as long 
as the difference Va = Wq;_l — ^ (a = 1, . . . , 25* -|- 1), the peak splits into as many peaks as distinct values of ujai^ 
exist. 

The study of the surface modes has also produced other results. We have found that the peak intensity of the 
surface phonon in the Laughlin liquid is proportional to l/m, where m is the inverse of the filling fraction. This 
was interpreted as if the bosons in the Laughlin liquid behave as "super-fermions" which, on average, occupy two 
states per particle. This is a particular example of the generalized exclusion statistics first discussed by Haldane [2^ . 
This feature of the spectral weight of the surface modes of the Laughlin state can be interpreted as an experimental 
manifestation of this exotic exclusion statistics. In addition, by analyzing the behavior of one-body density matrix 
near the boundary, we are able to obtain the occupancy of the orbitals with the highest angular momentum in the 
Laughlin state. For m > 1, we have argued that the behavior of the occupancy near the boundary forces a maximum 
in the occupancy before it decays to zero at I = Imax = m{N — 1), where N is the number of particles in the droplet. 
This maximum has its correspondence in the density profile, which should be a hallmark of the Laughlin liquid, and 
possibly of other vortex liquids as well. Last but not least, we have argued that the Laughlin liquid and the singlet 
states of vector bosons, should be very stable with respect to three-body losses. This is because their wave functions 
vanish when any two particles approach each other. The Moore-Read state that provides a good approximation to 
the state of total angular momentum L w N'^ /2 for scalar bosons, should also be very stable since its wave function 
vanishes when any three particles come to the same point. However, since this state is only a good approximation to 
the actual ground state, we expect it to be less stable than the Laughlin state, which is exact. 

In spite of the successes of the constructions presented in this work, there remain a number of important open 
problems. An important one is related to the accuracy of the description given here to describe the surface modes of 
vortex liquids of scalar bosons. The description seems to be well established for the Laughlin state, but for the other 
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stable states found at lower angular momentum it should be checked numerically. This can be done by analyzing the 
low- lying excitations of the states found in Ref . il^ We expect that the analysis performed here will encourage further 
studies in this direction, which may prove or disprove the conclusions of our analytical work. Another interesting 
issue that has not been addressed in this article is the possibility that not all the excitations predicted by the parton 
construction are indeed surface modes. This problem was numerically analyzed for fermion systems in Ref. >36i , where 
it was argued that the effective theory for the surface waves of states in the (fermionic) Jain sequence is a Wi+oo 
minimal model instead of the larger (in the sense that it has larger spectral degeneracies, but the same spectrum) 
m(1) X su{n)i model, which results from the parton construction used in the present work (see Ref. 36 and references 
therein). Similar comments apply to the other states of vector bosons which have not been considered in this paper, 
but which were recently constructed in Ref. To sum up, we can say that we have just begun to scratch the 

surface of the rich phenomena offered by this new setup for strong correlation and quantum Hall physics. 



The author wants to express his gratitude to Erich Mueller for a useful suggestion. He also gratefully acknowledges 
useful correspondence with Nicola Wilkin and helpful conversations with Ignacio Cirac, Marina Huerta, Alexander 
Nersesyan, Antonio Sarsa, Subodh Shcnoy, and Peter ZoUer. Financial support provided by the ESF programme 
"BEC 2000+" for a visit to the MPQ in Garching (Germany) is gratefully acknowledged. 

APPENDIX A: RELATIONSHIP BETWEEN PROPERTIES IN THE LABORATORY AND ROTATING 
FRAMES. THE DIPOLE MODE AND GENERALIZED KOHN'S THEOREM 

In EEC's of ultra-cold dilute gases time-of-flight and in situ imaging of shape oscillations, as well as Bragg scat- 
tering [s^, 113 ) have been used to measure the excitation spectrum. Light scattering in general, and Bragg scattering 
in particular, allow to access the dynamic structure factor 5(q, w) [s^ill^- Whatever the method of choice, here 
we shall imagine a situation where a time-dependent non axi-symmetric potential is added to the trapping potential 
in such a way that it drives the system slightly out of equilibrium by exciting surface modes. The response to the 
potential is measured and since we are interested in the linear response regime, we shall assume throughout that the 
perturbation is weak. All the measurements are performed in the laboratory (or lab, for short) reference frame. At 
this point it is convenient to recall the steps needed in the preparation of a rotating cloud. One of the commonly used 
methods consists in adding a time-dependent and non-cylindrically symmetric part, pretty much like the perturbation 
that we shall consider later, which imparts angular momentum to a conveniently cooled BEC. We shall assume that 
the cloud has been subjected to this "stirring potential" for some time, driven at some angular frequency Ud- When 
enough angular momentum has been deposited in the system, the external force is turned off and the could reverts to 
a cylindrically symmetric situation, where the Hamiltonian has the form discussed in Sect.^ Under these conditions, 
for an axially symmetric trap, the angular momentum along the axis, L, is a conserved quantity (i.e. it commutes 
with the Hamiltonian). This means that the total angular momentum imparted to the system will remain constant 
subsequently, rotating at a frequency fl ^ JlS^. Next we imagine that in this situation the system is perturbed 
again by switching one a weak time-dependent stirring force. In the lab frame, this perturbation is described by 
adding to the Hamiltonian the following term: 



In the previous expression, we have neglected any dependence on the axial coordinate. Furthermore, in order to excite 
predominantly surface modes, a multipole expansion of the external potential 5vLAB{'r,0,t) should not contain high 
multipolc terms. This means that it will be a smooth function of r near the boundary of the boson droplet, i.e. for 
r R. If we further assume that the external potential contains no monopolc mode, i.e. 
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where the last expression follows from the smoothness of the potential near the boundary of the droplet. The expression 
integrated over r is the boundary density operator, j{0). Hence 

ULABit)^ / de5vLAB{R.e,t)m. (A4) 

This corresponds to the perturbation discussed in the main text, but expressed in the lab frame. The rot and lab 
are related by the unitary transformation 

U{t) ^ e'■^"*/^ (A5) 

Notice that this transformation will affect only those correlation functions which are time-dependent (e.g. like the 
dynamic structure factor). Time- independent correlation functions, such like the one-body density matrix and its 
diagonal part, the equilibrium density profile, are unchanged when going from one reference frame to the other. To 
see how the dynamic structure factor computed in the rotating (rot, for short) frame is related to measurements in 
the LAB frame, which is where the measurement apparatus is actually placed, we need to work out the effects of the 
perturbation, Eq. ljA4|l . in lab. Starting from the time-dependent Schrodinger equation, 

indt\<^LAB{t)) = [HlAB + ULAB{t)] |$LAi3(0>- (A6) 

we find that, to linear order in ULAB{t), the wave function is modified to 

^ABit)) « |<fi^B(- y)) - J [\ dt' Ui^sit')\-l>iMt')), (A7) 

where the interaction representation has been used: 

\^\t)) = e'-f^^-*«*/'"'|$(t)), (A8) 

U\t) = eiHLABt/h^|^^y-~^Hr.ABt/h^ ^^g) 

Before going any further, it is convenient to recall that the Hamiltonians in the ROT and lab frames are related by 
the following expression 

Hrot = Hlab - nL, (AlO) 

where [L,Hlab] — 0, as the trap potential is cylindrically symmetric (this is not so after ULABit) is turned on). 
This expression is obtained by transforming the time-dependent Schrodinger equation, (|A6|) to the rot frame using 
Eq. I|A5|I . Indeed, it is the quantum analog of the classical relationship Erqt = E^ab — ^L. 

Assuming that, in the ROT frame, at time t = the system is in a state with probability given by the 

Boltzmann weight e~^" we compute the probability per unit time for large To that the system is found, also in 

the ROT frame, in a state at i = To/2 (we choose HroA^o) = E^^'^\<^a) and T|$q) = Lq|$q,), a = n,m. The 
calculation, however, is performed in the lab frame). The result can be written as 

= ^ E j ^\5vlab{R.I.^)? SROT{l.^-m- (All) 
Therefore, we can define the dynamic structure factor in the lab frame as 

SLAB{l00)^SROTil,UJ~m). (A12) 

where the dynamic structure factor, 



+00 /'27r 



00 



(/,^)= / dt / de'e^-'-^^'\j{e\t)m) = '^-^ye-^^"^'^\{m)W)^^ (A13) 



a, (3 



There is a simple explanation for the transformation law, Eq. ljA12|l . The shift w — > a; — ilHs precisely the frequency in 
the ROT frame of an excitation of angular momentum Til and energy fiuj in the lab frame. The other component in the 
dynamic structure factor are the oscillator strengths |(a|j(Z)|/3)p for transitions induced between different eigenstates 
of Hrot (and also of H^ab, since it commutes with L). Since the rot and lab frames are related by a unitary 
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transformation, which leaves these probabihties unchanged, only the delta function in lo is affected by the change of 
reference frame. 

Eq. (jAll|l implies that we can measure Slab{1,^) by cleverly choosing the external perturbation Svlab- It is 
important to point out that in the discussion in the main text, we have omitted any reference to the dependence of the 
intensity of the response on the droplet radius, i?, coming from SvLAsiR, l,'-^)- When analyzing the experimental data 
this fact must be taken into account since R, defined as the semiclassical rad ius o f the droplet, depends o n the effective 
filling fraction v. For instance, for a Laughlin droplet we found that R = y/roN i, and in general R = ^ N jv being 
the average density given by ly/irP. Thus, besides the factor v = p/{p + 1), with p a non-negative integer, discussed 
in previous sections, the strength of the peak in response to an /-polar potential Svlab ^ e*'^ is also proportional to 
i?'. This factor is especially important when analyzing the dipole mode (i.e. / = 1). This is because the energy and 
the peak intensit y o f this mode are given, for a many-body system confined in a harmonic trap, by the generalized 
Kohn's theorem |5lL Is^ l . This theorem follows from the fact in a harmonic trap the motion of the center of mass 
decouples from other motions of the system. Introducing the operators: 

N 

n = ^ [p, - Mm X r,] , (A14) 

J=l 
1 ^ 

i=i 

the Hamiltonian (in the ROT frame) splits into two commuting parts (cf. Eq. (0)), 

Hrot — HcM + Hr, (A-16) 

To study the excitation of the center-of-mass degrees of freedom it is convenient to use the set independent operators, 

At = Ik±^-z(c^^-|-r!)Ar(x + ir), (A18) 



M 

= +^{coJ^~m{X-^Y), (A19) 

ivl 



which have the fohowing properties: 



[Hrot,A^] = HluaA^ (A20) 

[Hrot,B^ = huBB\ (A21) 

[L,A^] = nA\ (A22) 

[L,B^] = -fiB\ (A23) 

where loa — lo±^ — il, i.e. the energy of / = 1 surface mode, which in the lab frame will have an energy equal to fiuj^. 
The theorem also predicts the existence of another Kohn mode, at lor — wj^-l-O in the ROT frame. Notice that although 
this mode has much higher energy in the ROT frame, it actually reduces the angular momentum by one quantum and 
therefore, application of the formula Erqt = Elab — ^L, yields an excitation energy in the LAB frame equal to uj±. 
Thus we conclude that the center-of-mass modes created by acting with the operators and B^ on the many-body 
ground state are degenerate in the LAB frame. However, since the system is in thermal equilibrium in the ROT frame, 
it is the energy in this system that enters in the Boltzmann factor, e~^" ^'^ . Hence, at temperatures T <C hfl, 
the probability of finding the system in the state ^^[0) is overwhelmingly larger than the probability of finding it in 

B^O). By their quantum numbers we can identify A^ ^ si — J2f=i ^ji ^'^d B'^ ^ di ^ SjLi ^ X^jLi %i ^'^d as 
we argued in Sect. IIIII the former yields a low-lying state whereas the latter does not. 
Let us finally compute the the peak intensity of the mode created by A^: 

f{l = l,T = 0) = \{l\N{X + zy)|0)|2 = ^|(l|(At - (A24) 



where |1) = |0>/^|[A, At] |. Thus we obtain 



f{l^l,T^Q)^^^Ne. (A25) 
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This is in agreement with the predictions of effective field theory once the dependence on R has been taken into account. 
For instance, consider a LaughUn droplet perturbed by the dipolar potential N{X + zy)e~*'^*, i.e. SvLAB{R,0,t) ~ 
i?e*^e-'"*. Using Eq. PlT|l and the results of Sect.Hlil 

fn = l^T = 0) = R^—^Nf (A26) 

TO 

where we have used R — V mN£. This result is identical to the one obtained by the generalized Kohn's theorem, as 
it should be. For states with several branches of surface phonons R^ = N£'^/i', and since the total spectral weight 
of Slab{1,^) for I = 1 equals v = p/{p + 1) {p integer), and taking into account that R^iy = Ni'^, we see that all 
the spectral weight must be concentrated in a single peak at uj — lu_i_. Therefore, all corrections to the energy of 
the surface modes must vanish for Z = 1, as discussed in Sect. HVl making all branches degenerate at this particular 
angular momentum. 



APPENDIX B: PLASMA FORMALISM AND VERTEX OPERATORS 



In this appendix we shall introduce the plasma formalism, which is a convenient tool to study some properties of 
quantum Hall wave functions. However, our main goal here is to establish a deep connection between this formalism 
and the correlation functions of a two-dimensional gaussian field theory [23, . As we show below for the Laughlin 
state, the TV-particle wave function can be obtained as the multipoint correlation function of N objects, known as 
vertex operators, plus another operator, which yields the gaussians omitted so far. The vertex operators are then 
related to the boundary field operators introduced in Sect. IIHI The construction shows explicitly the relationship 
between the low-energy description of surface modes, and the bulk wave functions. The results of this section are used 
in Sect.[3for the slightly more complicated case of vector bosons. Finally, we also show how to obtain the one-body 
density matrix, Eq. H2t)|) . 

Consider the wave function for Laughlin state 



N N 

|2 /ofl2 



. . . , zjv) - n - ^^o" n ■ (Bi) 

i<j — l i—1 

The (unnormalized) probability density can be written as 

P™(ri,...,rAr) = |$„(zi,...,zjv)P = exp[-/3C/(ri,...,rjv)], (B2) 



where /3 = 2to and 



^ TV N 
i—l — l 

is the potential energy of a one-component classical plasma in two dimensions. The plasma consists of N identical 
particles with unit charge moving on a uniform neutralizing background of charge density po = I/tott^^. The plasma 
is neutral if the background takes up a disk of radius R — ^ {Trpo)^^N — \/mN£. To understand the analogy, let us 
recall that in two dimensions the Coulomb potential due to a unit charge reads 

g{r — r') = — log |r — r'| + const. = — log \z — z'\ + const. (B4) 

In addition, a uniform background of charge po creates a potential at r (with |r| < R) given by 

ffBG(r)=/ dr'g{r-r')po = -^\z\^ (B5) 

as can be easily checked using Poisson's equation: V^(7BG(r) — — 27rpo- Putting everything together, Eq. (jB3|) is 
obtained. 

It will be useful to rewrite (|B2|I as follows 

P^[p]-^^y (B6) 
gauss |_^J 
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where ^gauss[p] is the functional integral 

^gaus.[p] = / [du{r)] ^-PJdr[^(Vuir)f-^pirMr)]^ ^g^^ 



where p(r) — — X^i^i '^(i" ~ ^i) + Po'&{R — |r|) is the total charge density of the plasma {d{r) is the step function). The 



above expression can be shown to be equivalent to Eq. (|B2(I by shifting 

u{r) = u'{r) +i dr' g{v - rXr'), (B8) 



which cancels the term linear in u{v) in ljB7p and yields the following expression for the potential energy of the plasma 

C/(ri,...,r^) = C/[p] = i j dvp{v)g{v-v')p{v'). (B9) 



Substituting p(r) defined as above, this expression reduces to ljB3p (the charge self-energies can be removed by 
redefining g{v) = — log [(|rp + a^)/a^] and carefully taking a ^ at the end). 

We have taken so much work to show that P(ri, . . . , rjv) can be written as the functional integral (|B7() because this 
form allows us to make a very useful connection. Indeed Eq. (|B7|1 is the "generating functional" of a two-dimensional 
field theory: For arbitrary p(r), it generates all correlation functions of a gaussian field theory, which describes, e.g., 
the classical statistical mechanics of a two-dimensional superfluid or the quantum mechanics in imaginary time of a 
one-dimensional quantum fluid I3- K we specialize to p{v) — po'&{R — |r|) — Yli'i=i ^{^ ~ ^i)^ then Eq. HB7|I generates 
A^-point correlation functions of the form 

P(ri,...,rAr) = (e-2'"""(''i)...e-2™™('""'e i-i<« ), (BIO) 

If instead of regarding this correlation function as a statistical average in two dimensions, we think of it as a quantum 
expectation value in one dimension, the objects e2mm(r) j-gpj-gggni; operators, termed "vertex" operators. The field 
operator u{r) obeys the equation of motion V^u(r) — 0, i.e. the Laplace equation. This also corresponds to the 
equation of motion of a phonon field in imaginary time. Recalling that any solution of Laplace's equation can be 
written as the linear combination of an analytic and an anti-analytic part, let us write 

u(r) = [^(z) - ^(z)] , (Bll) 

where z = x + iy and z = x — iy and the normalization \/ ^/m has been introduced for later convenience. Next, notice 
that the above field operators obey the (first order) equations of motion: 

dz(p{z) = i ((9a; - idy) q)(x - iy) = 0, (B12) 

d-Mz) = \id^ + ^dy) <fi{x + ly) = 0. (B13) 

The structure of these equations resembles very much the equation of motion for the chiral phonon field 4'{9) intro- 
duced in Sect. IIIII Recall that there it was found that (f){0,t) — (f){9 — vt), which in imaginary time t = it obeys 
(^v~^dr + ido) <j){vT — i9) — 0. In an analogous way as we did for (j){9) in terms of e'^, we can expand the fields (/? and 
Cp 'm a series of z and z, respectively. These expansions read: 

ip{z) = _<^„ _ log z -f ^ i= [z-'c{l) + z'ct (0] , (B14) 



/>0 



^{z) = ^o + iglogz + ^i= [z-'c(0 + z'c^(0] , (B15) 



i>0 



where [g, (po\ — i, [c{l), c^(^')] — ^i-i' , and similar expressions for g, ipo and c(Z), c^(/). Returning to the vertex operators, 
which are our main interest in this appendix, we notice that, in operator formalism, it is necessary to normal- order 
them as follows: 

Vg(z) = : e*'3'^(^) : = e'"'^^'' e'^'Ji"^^ e"^^'>o e"'S'>o ^''^^^^^ (B16) 
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and a similar expression holds for the vertex operators of ifiz). This is needed in order to avoid divergences, which 
are related to the "charge self-energy" problem encountered above when using the functional integral approach. With 
these provisos, the probability distribution of the Laughlin state can be written as: 

P„(ri, . . . ,rAr) = . . . ,^jv)$„(zi, . . .,zn), (B17) 

where 

-\-i^/m J d^Zip(z)pa 

^^{z,,...,zn) ^ {V_^{zi)...V_^X^n) -.e i^K- :), (B18) 

— ^^/7n J d? z (p{z)pa 

<i>ra{zi,...,ZN) ^ {V+^{z,)...V+^{zn) -.e i^i<« :>. (B19) 

Therefore, the probability distribution for the Laughlin state splits into the product of two correlation functions of 
the gaussian theory. Let us now demonstrate that each of these equals either the Laughlin state or its anti-analytic 
version (i.e. its complex conjugate). To this end, it is useful to further investigate some of the properties of the vertex 
operators. For example, their two-point correlation function reads: 

{V,iz)V^,iz'))^{z-zr''\ (B20) 

as it can be shown from the identity e^e^ = 6*^6^6"^ for any two operators with a c-number commutator C = [A, B]. 
The generalization to multi-point correlation functions can be obtained in a similar fashion: 

{Vp, {z,)... Vp, {zn)) = n - 'i)^'^' 
i<j 

provided that X^i^i Pi — ^- Using this result, it can be shown that 

N 

$„(zi,...,z^) = e™(^--'^~) X{{z,-z^r\le-\'-\"'^'\ (B22) 

i<j 2—1 

i.e. the Laughlin state. The phase factor e™(^i' - '^™) cancels out after taking the product with the anti-analytic part 
p{zi, . . .), and we shall not consider it any further. 

To end this discussion, let us make firmer the connection to the low-energy description of surface modes, which we 
have loosely established above. This can be achieved by a re-examination of the field operator found in Sect. lIIII There 
we showed that the boson field operator is related to the operator e='=""'^(*), which can be shown to be proportional 
to a vertex operator. Indeed, if we normal order it properly, we obtain (a is a short distance cut-off): 

^±im<j)(e) _ ^m/2g-ime/2 . ^±im(j)(e) . ^ (B23) 

We next notice that by comparing the mode expansions, Eq (|B15() and Eq. H21|l. one can identify 0(0) = (f{z — 
e~'^^)/^/m. Thus we arrive at 

■^{9) = Ae+'"'(^-^'>'^V+^{z^e-''^), (B24) 
¥{9) = Ae'"'^^V_^,{z^e"'')., (B25) 
from which one can readily show that the one-body density matrix at the boundary, 



G{0 - e') = (*^(6')*(e'')) = const X e-«'"'°(«-«') 



(B26) 



as follows from Eqs. (|B24IB25p . and IIB20|) . Therefore, if we define the boson operator in the LLL as '^^q{z) = V_^/j^{z), 
the Laughlin wave function can be written in the very appealing way. 

This is the main result of this appendix. 



I>„,(zi,...,zjv) = (*|j(zi)...*|3(zAr) :e i--r<« :). (B27) 



[1] A discussion of the role of broken Galilean and f/(l) Gauge invariance in classical and quantum fluids can be found in 
A. M. J. Schakel, Boulevard of Broken Symmetries, cond-mat/9805152 As explained by Schakel, another consequence of 
broken Galilean invariant is the existence of sound waves, which are the Goldstone modes. 
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